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ABSTRACT 


The relationship between beam coupling losses in multiple beam antennas and 
pattern parameters such as beamwidth and sidelobe level is investigated. A computer 
algorithm to calculate the coupling coefficients for arbitrary number of beams is 
implemented using MATLAB. Gain-loss data for various array amplitude distributions 


is presented. 
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I. INTRODUCTION 

An important component of all communications, radar, and electronic warfare 
systems utilizing the radio-frequency portion of the electromagnetic spectrum is the 
antenna. The antenna is a transition structure that serves to minimize the signal loss 
between a guiding structure and the transmission medium. It must efficiently radiate 
power when transmitting and efficiently collect power when receiving, and is usually 
called upon to give directional information as well. With so much reliance on the 
electromagnetic spectrum for detection, targeting, and communications, it 1s imperative 
that technical personnel understand the roles and design constraints of antennas. 


There are three principal types of antenna systems in use: 


1. Reflectors 
2. Lenses 


3. Arrays 


Although multiple beams can be formed using any of these antennas, the analysis 
presented here will concentrate on arrays. An array is simply a collection of small 
antennas, generating a directional beam using constructive and destructive interference 
[Ref. l:p. 204]. Array antennas are widely used in_ radar, and 
communications applications. Naval applications include the AN/SPY-1 radar for the 


AEGIS system, the AN/AWG-9 F-14 radar, and the receive multiple beam antennas for 


the AN/SLQ-32 electronic warfare system. Additionally, a conventional reflector 
monopulse radar consists of a 4 feed horn assembly arranged in a planar array. 


Arrays are particularly well suited to the radar role for the following reasons: 


1. High data rates are possible with electronic beam scanning. 
2. Multiple beams can be generated. 


3. Lower sidelobe levels are possible for Low Probability of Intercept Radar 
(LPIR). 


4. Solid State technology can be efficiently utilized. 


Before array technology, increasing directivity meant increasing the electrical size of the 
antenna relative to the wavelength of the transmitted signal. The problem of increasing 
the electrical size of an antenna was particularly exacerbated at lower featenc es Arrays 
were necessitated since it was not feasible to generate a single element with the required 
dimensions. For example, at a frequency of 30 Megahertz (MHz), the corresponding 
wavelength is 10 meters (= 30 ft). A two wavelength antenna would be 60 feet long. 
Alternatively, the same directivity could be obtained using 2 or 3 quarter wavelength 
antennas each 2.5 meters long as an array, with half-wavelength spacing. The efforts of 
pioneers like Schelkunoff led to new advances in antenna design using array techniques 
and polynomial methods. 

In addition to electronically scanning the array with a single beam, it may be 
advantageous to receive multiple beams simultaneously. Multiple beam antennas offer the 


capability to scan sectors rapidly with little reduction in antenna gain. As mentioned, 


many systems currently in use by the United States Navy take advantage of multiple 
beam antenna technology, including the AN/SPY-1 radar, the AN/SPS-48 radar, and the 
AN/SLQ-32 Electronic Warfare suite. 

The desire to optimize multiple beam systems has led to new synthesis techniques. 
Advances in computing ability and efficiency have led to various methods for generating 
multiple receive beams. These beam forming networks, with a matrix interpretation are 
well suited for implementation using array elements as 
feeds.[Ref. 2:p. 241] The advantages offered by phased array radar 
systems (fast scanning, low sidelobes, ECCM capability) combine with the advantages 
of multiple beam systems to further enhance overall system performance. Currently, most 
modern systems that employ multiple beam antennas rely on some type of array to 


implement the design. 


A. MOTIVATION 

High performance radar systems need narrow “agile” beams to accurately track 
targets. Many simultaneous beams are desired to track multiple targets. The necessity for 
high data rates common to most current weapons control systems makes electronic beam 
scanning necessary (as opposed to mechanical beam scanning). Given that premise, 
phased array antennas are the antennas of choice in most of these new systems. 

Recent advances in Monolithic Microwave Integrated Circuits (MMIC) technology 
have expanded the number of RF components available in a small package. Currently, 


emphasis has been placed on integrating all of the RF devices necessary for a system (i.e. 


phase shifters, radiating elements, low noise amplifiers, etc.) onto a single substrate. 
These technological advances have produced compact arrays with a large number of 


multiple beams. 


B. PROBLEM STATEMENT 

Intuitively, we would expect interaction between the various beams in a multiple 
beam system. For example, if one of the beam ports is excited, some of the energy will 
couple to the remaining beam ports and appear as a received signal. This coupled energy 
results in a loss in gain that would not be present if the antenna only had a single beam, 
and occurs even if there is perfect isolation between the beam feeds. This interaction has 
been investigated by many researchers, most notably Seymour 
Stein.[Ref. 3:p. 548] These beam coupling losses are | an important 
consideration in the design and analysis of multiple beam antennas. There is a tradeoff 
between these losses and the beam crossover level. High crossover (i.e. closely spaced 
beams) are particularly desirable from a systems perspective. However, the beam 
coupling losses at these crossover levels may be excessive. In receive antennas, an 
adequate signal-to-noise ratio is maintained using low noise amplifiers (LNA’s). The 
location of the amplifiers in the antenna feed and their required gain are influenced by 
the beamforming losses. Thus prediction of the beam coupling losses is essential as the 
number of multiple beams, and the crossover level increases. Figure 1 shows the 


radiation pattern for a typical multiple beam antenna array using a Taylor distribution 





with 25 dB sidelobes and n=5. In this figure, there are 3 beams and the beams 


crossover at the 3 dB point. 


C. APPROACH 

Typically, beam coupling losses are estimated from curves such as those by Stein 
or from crude estimates obtained by considering only two beams and extrapolating to a 
higher number of beams. These techniques are effective for a few beams, but accuracy 
rapidly drops as the number of beams increases. Once an antenna development model is 
available, beam coupling losses can be measured, but this technique is costly and time 
consuming. 

Stein has formulated the loss calculation in terms of a cross correlation integral and 
antenna scattering parameters. His formulation has reduced the problem to solving a 
matrix eigenvalue equation. Flexible software computing packages such as MATLAB” 
or MATHEMATICA’ can be used to solve the equations for a wide range of geometry 
and excitation with a single algorithm. 

The goal of this thesis is to compute the beam coupling losses in a multiple beam 
linear array as a function of the number of array elements, the element spacing, the beam 
crossover level, and the sidelobe level. For computational convenience it is assumed that 


the excitation coefficients for all beams are the same (the case of identical beams). 


D. CHAPTER OVERVIEW 

A breakdown of contents by chapter is provided here. 

Chapter II lays the requisite mathematical foundation in antenna theory and gives 
a derivation of the fundamental relationships in array design. 

Chapter II gives the mathematical foundation of the theory of beam coupling losses 
using a scattering matrix representation. The adaptation of Stein’s work to linear arrays 
is also presented. 

Chapter IV presents the summary of the calculated data for the two and four beam 
cases, and looks at the relationship between beam coupling losses for a given geometry 
and the number of beams. 

Chapter V summarizes the work performed and gives the conclusions. 

Analytical details, raw data, and computer source codes are included in the 


appendices. 


Il. ARRAY FUNDAMENTALS 

One definition of an array is "... an assembly of radiating elements in an electrical 
and geometrical configuration. "[Ref. 1:p. 204] Arrays produce radiation patterns using 
the principles of constructive and destructive interference. Destructive interference is 
used to suppress radiation in unwanted directions. The interference patterns are a function 
of the weighting of the individual elements. The weighting can in general be complex, 
and determines the beamwidth, sidelobe level, and scan _ direction 
[Refs4:pe i: The design parameters that influence the characteristics of 


the radiation pattern are: 


1. geometrical configuration, 
2. relative element separation, 
3. element amplitude weight, 
4. element phase weight, 


5. radiating element pattern [Ref. 1:pp.204-205]. 


Maxwell’s equations are used to derive the expression for the radiation pattern. 
Two assumptions that apply to antenna problems are 
1. All quantities have a time-harmonic dependence. Phasor notation will be used 
with the e’ suppressed. 


2. Antenna surfaces are perfect electrical conductors (i.e., only J, ). 


The electric field at a field point P(x,y,z) due to the current Hoey ,Z’) can be expressed 


as [Ref. 1:pp.86-91] 





E(x,y,2) = -joA - WV - A), (1) 
| Ho &o 
where 
A(x,y,2) = Bec oy ds’ 
(x,y) real [F,@'y’2)— ae 
je a + (y-y/? + — 


S denotes all of the antenna surfaces; for an array it becomes a collection of identical 
discrete surfaces. Combining (1) and (2) and simplifying leads to the expression for the 


electric field 








x ~jkr (3) 
E(xy,2)=-> ale (xyz) —— v vl effie's's asp 


In the far field Ro, and the second term in (3) can be neglected since it will involve 
terms decreasing as 1/R’ and 1/R’. Using the geometry of Figure 2, and making the far- 
field approximation that R and r are parallel gives |r| =~ |R| -R’ - R. The electric field 
at P becomes 


EC eile if if J, - VR) Re e FRR ds (4) 


where R is a unit vector in the direction of point P. 
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Figure 2. E Field Calculation Geometry (Far-Field) 


Now that the expression for the electric field has been derived, it 1s necessary to 
extend this result to multiple elements. Consider an array consisting of N identical 
surfaces, and define a reference element with surface S, as shown in Figure 3. The 


expression for the electric field of the array can be written as 


E(x,y,2) a Zotoe wey. “il Js - (J, RR] elas (5) 
0 


where c; is a complex constant that relates the current on the i” element (ia) to the 
current on the reference element (J ,). By factoring the complex exponential term in the 
integral, (5) can be rewritten as 


Fee [I Pos ee ta 


0 


5 a cE: 
ECs) a)a) eae. 
i=] 











ea) ce) 


From (6), it can be seen that the total radiation pattern for the array is the product of the 
element radiation pattern (EF) and the array factor (AF, the term involving the 
summation). This is known as the principle of Pattern Multiplication. 

For radar applications, the antenna gain, beamwidth, and sidelobe level are all 
important design parameters. High gain increases radar range and a narrow beam 
improves target resolution. However, low sidelobe levels have currently become a 
priority for clutter reduction and electronic counter counter-measures (ECCM). 
Unfortunately, low sidelobes are achieved at the expense of both gain and beamwidth. 


Thus, there is a relationship between beamwidth, sidelobe level, and gain that forms a 
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Figure 3. Array Element Geometry 





classic engineering tradeoff. The sidelobe level is set by the proper choice of the 
magnitude of the constants, c;. Table 1 gives a comparison between the beamwidth and 
sidelobe level of several commonly used amplitude windows. 


TABLE 1. BEAMWIDTH - SIDELOBE COMPARISON 


Window (|c;| =a) 


a 


as ——_—— 


Sidelobe Level (dB) HPBW (VL) 


; * 
Taylor (20 dB n=4) : ize 


Hamming | -42. 1.302 


The Taylor weighting is used extensively in radar applications, and is specified by two 
parameters: sidelobe level (SLL) and n. For a given sidelobe level, it offers minimum 
beamwidth. Once either gain or SLL is specified, the other quantity is established. [Ref. 
1:p.680] 

Thus far only the effects of amplitude weighting the array elements have been 
discussed. However, as stated previously, the element weighting can be complex. The 
effects of phase weighting will now be addressed. Let the N isotropic radiators lie along 
the z-axis with element spacing d as shown in Figure 4. According to (6), the array 
factor is the linear superposition of the contributions of each individual element. In terms 


of the quantities in Figure 4, the array factor is 





Figure 4. N-Element Array Far-Field Geometry 


j(kdcos® + B)(n- X*1) 


N 
AF(8,6) = ¥> a,e a, (7) 
n=] 


Implicit in equation (7) is the representation of the complex weights as 
cmm= ac ibm (8) 


where a, 1S the amplitude weight and 8 is a progressive phase shift between array 
elements. If a,=1, the array is a uniform array, and the sum becomes a geometric series 


that can be reduced to the closed form 


sof 
_i¢_/| (9) 
2) 


“(3 


where y=kd cos 6+8. If 6=0, the direction of maximum radiation is perpendicular to 


1 
AF(p) = W 
the axis of the array elements, and the beam is referred to as a broadside beam. The 
pattern for a 10 element uniform array with d= )/2 is shown in Figure 5. If 6 =kd, the 
direction of maximum radiation is in the direction of the array axis, and the beam is 
referred to as an endfire beam [Ref. 1:p. 218]. 

The advent of ferrite materials and high speed switching circuits have made 
scanning radar arrays possible. By applying a progressive phase shift to equation 1, the 
far field radiation pattern is shifted by an amount determined by the linear phase 
variation across the array. The beam peak location @) for the amount of phase shift 6 is 


determined from 
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Figure 5. 10 Element Array Factor - Broadside Scan 
(15°/sector) 
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Bp = -kdcos6, . (10) 


Inserting this result into (7) and simplifying leads to 


N » 
AF(6) =) Ic, (a aaa = cee (1 1) 
n=] 


where 6, is the desired scan angle. Broadside and endfire correspond to 6,>=90 and 6).=0 
or 180 degrees, respectively. For radiation in a particular direction between 0 and 90 
degrees, the appropriate phase 6 must be selected according to equation (10). The array 
can then be made to scan in any desired direction, hence the term "phased array”. [Ref. 
iep.2 20] 

An example of a radiation pattern for a 10 element uniform array scanned to 45 
degrees is shown in Figure 6. For most radar applications, the scan is predominantly + 
60° from broadside. Furthermore, only radiation into one hemisphere is desired. To 
achieve this, the linear array is placed above a ground plane. If the ground plane spacing 
is properly chosen, the net effect 1s to leave the pattern above the ground plane 
essentially unchanged and eliminate the radiation in the back hemisphere. Linear arrays 
above a ground plane can be modelled as a planar array radiating in free space by the 
method of images.{Ref.1:p. 137] The pattern for a 10 element uniform array over a 
ground plane scanned to 45 degrees is shown in Figure 7. The gain will increase by a 
factor of approximately 2 if the power into the antenna is kept the same, since all of the 


power is confined to a single hemisphere. Current systems that employ phased array 
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Figure 6. 10 Element Array Factor - 45 Degree Scan 
(No Ground Plane) 
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Figure 7. 10 Element Array Over Ground Plane - 45° Scan 
(15°/Sector) 
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technology include the AEGIS AN/SPY-1 series of Naval Radars, the PATRIOT Air 
Defense Weapons System, and the PAVE PAWS Early Warning Radar System. 

The final array design parameter of interest is the interelement spacing, d. When 
scanning an array, multiple directions of maximum radiation within the visible region 
(0<6<180 degrees) may arise. These additional maxima are not desirable and are 
referred to as grating lobes [Ref.5:p. 114]. The choice of d_ influences 
the number and location of these grating lobes. The grating lobes appear when y is an 


even multiple of 27. Accordingly, the locations of the grating lobes are given by 
cos 8, = cos 8, + m (12) 
m =+1,+2,-- 


To ensure that no grating lobes appear in the visible region requires that cos 6, > 1, and 


thus the limitation on spacing becomes 


A | (13) 


1 + |cos6, | 


When 6, = 0, cos @ =1 and when 6) = x, |cos 61 = 1. Therefore if d < X/2, no 


grating lobes will appear for any scan angle @) [Ref. 6:p.144]. 
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HI. MATHEMATICAL FORMULATION OF BEAM COUPLING 


A. INTRODUCTION 

An antenna can be modelled as a multiport junction with scattering parameters S;. 
The S; are computed from a knowledge of the beam radiation patterns. Once the beam 
radiation patterns and the S; are known, the radiation efficiency can be determined by 
solving a matrix eigenvalue equation. This solution technique will be specialized to the 
case of a phased array antenna system. The model is based on the work of Stein [Ref. 


2], and his notation will be used. 


B. SCATTERING MATRIX FOR N-PORT NETWORKS 

For electrical networks at low frequencies (f < 20 khz), the standard lumped 
parameter circuit models work well in describing the behavior of the circuit. The circuit 
can be described by voltages across and currents through each element. The voltages and 
currents are related by the element impedances and admittances. However, relationships 
for low frequencies are not valid for microwave frequencies, since the phase of the 
electromagnetic wave may change significantly over the length of the 
element.[Ref. 7:p.1] The lack of the standard circuit models at microwave 
frequencies makes measurement of voltage and current extremely difficult. For this 
reason, at microwave frequencies it makes more sense to work with the incident and 


reflected fields in a network. The use of incident and reflected (or scattered) fields leads 
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to an equivalent network formulation, called the scattering matrix representation. [Ref. 
7:p.221] Scattering matrices mathematically describe the behavior of waves at each port 


of a device. For an N-port network, 


y2 So, Sy. Soy [2 (14) 


wn} [Sxz Sw2 Sew} |) 


where y denotes an outgoing wave and x an incoming wave. The elements S;; are called 


the scattering parameters of the device. In matrix form 
y=Sx . (15) 


The S, are found by driving port j with an incident wave, and measuring the outgoing 
wave at port 1, with all other ports terminated in matched loads. 

It is possible to visualize a multiple beam antenna as an N-port network consisting 
of N input feeds and N output beams. If the network is reciprocal (as is usually the case 
for most antennas that transmit and receive), the scattering matrix will be symmetric. 
[Ref. 7:p.224] The scattering matrix can be related to other matrix formulations, such 
as the impedance or admittance matrix. As with the latter two, the scattering matrix 
provides a complete representation of antenna performance[Ref. 7:p.223]. Figure 8 gives 
a conceptual picture of a multiple beam antenna system represented as an N-Port 


network. 
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Figure 8. N-Port Network (after Stein) 
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C. BEAM COUPLING FACTORS 

Using the notation in Stein’s paper [Ref. 2], and treating the antenna assembly as 
an N-port network, the incident fields at each port are represented by the vector x, and 
the reflected fields at each port are represented by the vector y. The x, and the y, are 


related via the scattering matrix 


N 
Ba aeer » Sim*m (16) 


The incident waves give rise to far-field beam patterns and the far-field radiation patterns 
are normalized using an autocorrelation integral over the spatial angular domain as 


follows 
i [rR (8,6)-R,(0,6)sinO6d6db = 2, , (17) 


where R,(6,@) is the beam pattern of the k“ beam ,y, is the impedance of free space 


(1202 (), and the “ denotes complex conjugate. The electric fields are of the form 


— ak kr 
E (0,0) = a, R,(0,6) — | (18) 


The term q, iS equivalent to the square root of the radiation efficiency and r is the 


distance from the local antenna coordinate system origin to the observation point. 
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The beam coupling coefficients are given by 


] 2n fT =. = ‘ 
B,, = ale I R, (8,6)-R(8,) sin6 dodo . (19) 


The coupling coefficients possess symmetry in k and j (6,;=6,; ) and are a measure of 
the beam overlap. Since the R,(6,¢) are vectors with complex components (i.e, each 
having a magnitude and phase), it 1s not just the area common to the polar diagrams of 
the multiple beams. The complete overlap of beams k and j implies 6,,=1; if 6,;=0, the 
beams are orthogonal. For uniform illumination, orthogonality occurs when beams cross 
at approximately the 4 dB points. 

Having established a measure of coupling, it is natural to investigate the 
implications of conservation of energy. If the antenna system is treated as a linear 


network, superposition can be used to compute the total radiated electric field. 





; ak a 20 
E(0,9)= | D7 x9, R(86)| — us) 
k=1 
Using equation (16), the incident total power P.,. is 
Pye, x 
me » —- (21) 


=e 
where the * indicates a conjugate transpose. Similarly the reflected power is 


P.=x'S'tSx . (22) 


ref 


The total radiated power, P,,,, is defined as 


jhe, 


P_, = 5 [77 no E'@.6) E(,6)sinddbdo (23) 


Appropriate substitution and simplification leads to an alternative expression for P.,4, 


N WN 
Prot = Dy dy Fe Ue By Fj; 


k=1 j=1 (24) 
N oN 
=), wars 
k=1 j=l 
Conservation of energy now requires that 
gene? Sr et oe (25) 


where equality holds only if the network is purely lossless. 


D. SPECIALIZATION TO AN ARRAY 


1, COUPLING COEFFICIENTS 


The adaptation of Stein’s theory to linear arrays as presented here follows 


original work documented in reference [8] [Ref. 8:p. 1]. The first step 
to determine the normalization constant for the k” beam in the array. The array element 
spacing is d, and the amplitude distribution is given by the coefficients {a,}. The 


magnitude of the normalization constant for the p" beam is defined as 


age | kdcos 6, ( I-5 a 
jkdcos® (m-n) 92 
> > atime ’ z ( 6) 





If m=n, L,,=2kd; otherwise, 
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I ee a | | (27) 
Suit (n - m)kd 

The beam coupling coefficients are computed using (19). Using the notation 
in Chapter 2, the array factor is substituted for the beam pattern, R,. Neglecting the 
element factor will not have a significant effect if the radiating elements are well matched 
over the range of beam scan angles. Thus, 


=f" es P AeA Sa ane P= 0 Ltud (28) 
1 


a= 


with Q(m)={2m-(N+1)}/2. The a,, and a, are the amplitude coefficients for the beams, 
the u, and u, are k times the phase shifts to scan the p” and q” beams, and u is the z 
direction cosine times k. Performing the integration and simplifying leads to the final 


expression for the coupling coefficient 


- N+] N WN 
J 7) (up -u9) j(u,m-u a 
1 € aa, é es 
B = m=1 n=l (29) 
Pq N N . AN ON ae 1 
Ju,(m-n 2 Jus(m-n 2 
Ey aero aay 
m=l1 n=l m=] n=] 








2. RADIATION EFFICIENCY AND COUPLING LOSS 
The calculation of the beam coupling loss is based on conservation of energy. 
Stein has shown that the problem can be cast into the form of a matrix eigenvalue 
equation. From the results of section (C(1)), the conservation of energy equation can be 


written as 


pa 


xix > x" SUS coe (30) 


Both of the matrices, T and S'S, are Hermitian and positive semidefinite. Thus, from 
linear algebra, each has N positive real or zero eigenvalues. Each matrix also has a full 
set of linearly independent eigenvectors. It is possible to construct a matrix, U, from the 
normalized eigenvectors described above that will diagonalize I and S'S. If U is used to 
diagonalize I’, the transformation is a similarity transformation and has the mathematical 


form of 


UTU=y , (31) 
where y is a diagonal matrix. U is also a unitary matrix, since it has the mathematical 
property 

UU =I , (32) 
where I is the identity matrix. 
Returning to the energy conservation equation, (30) can be transformed using 
the matnx U. Let x=Uz; substitution of this result into equation (30) gives 
zizg > zt'UTS'SUz + z'U'TUz . (33) 


The transformation from x to z is equivalent to exciting the antenna with a new set of 


voltages. Since x was arbitrary to start with, the new vector z represents an equally 
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arbitrary set of voltages. Rewriting (33) gives 


zid-y)z2z'(U's'sUWz , (34) 
in which U'S‘SU is also Hermitian and positive semidefinite. 
The characteristic equation of (34) is 


U- yf =0 | S>) 


Most standard mathematical software packages contain routines that can efficiently solve 
this eigensystem. The result is N eigenvalues (A,,A3,--,Ay ) that are related to the beam 


efficiencies ( lal’, lao |’, wat ln]? ) by 


] 
Ye = Wi Bua = — (36) 
At 
Because 6,,=1 by choice of normalization, 
] 
np es | 37 
Aj 


Since |q,|° < 1 from conservation of energy arguments, the upper bound on the 


radiation efficiency is given by 


] 
ie = (38) 
oe 
For the special case of identical beams, | q, | =|q,| =--=|qy|, the determinant 


for the radiation efficiency factor becomes 
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(1-A) Bis . Bin 


B,, (-A) - Boy re (39) 


Bay B vy. as (1-A) 


Using the symmetry of the coupling coefficients, allows the equation to be wnitten as 


(1-A) By By 


B., (1-A) aa Dae = (40) 


| Bayz B ye ae (1 -A) 
The evaluation of the determinant analytically proceeds by cofactor expansion. Expanding 


and collecting terms leads to the final result of 


INA AN +A AN? 4 4Ay AtAy = 0. (41) 


The largest root of this equation, X,,,,, 18 then substituted into (38) to obtain |q,,,,|°. The 
detailed solution for four beams is given in Appendix B. The results from Appendix B 
were used to verify the computer code which is based on the more general eigenvalue 
solution of equations (35) through (38). Note that identical beams can only be achieved 
where there is perfect symmetry, such as a circular array or a spherical array. However, 
the beams away from the edges of a large cluster of beams can be considered 


"approximately" identical. 
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A. 


IV. DATA SUMMARY 


INTRODUCTION 


The beam coupling loss data was computed for various amplitude weights applied 


to the array to control beamwidth (hence crossover level) and sidelobes. The amplitude 


weights used were 


B. 


ee 


5. 


6. 


Uniform 

Cosine on a 10 GB pedestal 
Cosine on a 15 dB pedestal 
Taylor (20 dB, n = 4) 
Taylor (30 dB, n = 5) 


Taylor (40 dB, n = 6) 


These distributions were chosen because of their practical use in radar and 


communication systems, and because of their ease of computer implementation. 


MATLAB has a built- in routine to calculate uniform weights, the Cosine and Taylor 


distributions were generated using author-created functions. Figure 9 shows samples of 


the various amplitude windows. 
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Legend 





—— Uniform 
—- Cosine (10 dB) 
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Figure 9. Amplitude Weightings 
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B. 2 BEAM DATA 

This section presents the summary of data for the 2-beam array scenarios. The 
complete raw data is presented in a set of tables located in Appendix A. 

The data summary consists of parametric plots of the coupling loss. With two 
simultaneous beams, the coupling loss is related to the crossover level through the 


coupling coefficient (8) between the beams 





q? = | (42) 


A plot of coupling loss vs crossover level for the 20 element array is shown in Figure 
10. The plot contains the data for all six amplitude weightings. Figure 11 shows a similar 


plot for a 100 element array. 


C. 4BEAM DATA 

Following is a summary of data for the 4 beam array scenarios. The complete data 
set 1S presented in tables located in Appendix A. As in the 2 beam case, the data 
summary consists of plots of the coupling loss. With four simultaneous beams, the 


coupling loss is related to the crossover level by 


2 I 
—, (43) 
A 


max 
where A,,,, 18 the largest eigenvalue of the 8 matrix. Plots of radiation efficiency vs 


crossover level for the 4 beam, 20 element array are shown in Figure 12. The plots 
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Figure 10. 2 Beam Coupling Loss - 20 Element Array 
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Figure 11. 2 Beam Coupling Loss - 100 Element Array 
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Figure 12. 4 Beam Coupling Loss - 20 Element Array 
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contain the data for all six amplitude weightings, and the corresponding plot for a 100 


element array is shown in Figure 13. 


D. MULTIPLE BEAMS AND GAIN LOSS 
One important result of multiple beam antenna research is the investigation of gain 
loss due to beam coupling. According to Hansen, the gain for a linear array is 
(Ref. 9:p. 20}. 
N Z 
zB 


@ ES ere (44) 
N 


N 
> nGm sinc (n-m) kd 
n=] m=] 


If the interelement spacing, d, is \/2, equation (44) reduces to 


a 3 = (45) 
a, 


which can also be expressed as 


| 


(46) 
/ 


G=Nn 


Bs 
where 7, 1s the aperture efficiency (1.e., gain loss due to non-uniform amplitude weights). 
This is the gain of an isolated lossless antenna with excitation coefficients {a,}. (In 
standard antenna terminology, equation (46) gives the directivity.) In practice, other 
losses will further reduce the antenna efficiency. The losses include transmission line 


losses, mismatches in the feed, and the coupling losses if multiple beams are formed. 
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Figure 13. 4 Beam Coupling Loss - 100 Element Array 
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The overall antenna gain is determined by the product of the individual efficiencies. When 


only tapering and beam coupling losses are considered for beam k 


G=Nnlal (47) 


As an example, the aperture efficiency of a single beam array with a 30 dB Taylor 
distribution is .8733. When a second beam is added with a 3 dB crossover level, there 


is an additional .7258 reduction in gain. For N=20, 


G = (20) (.8733) (.7258) = 12.67 (48) 


11 dB 


Thus the beam coupling loss is 1.39 dB. It is clear that any design incorporating multiple 
beam techniques must account for the coupling losses to assure proper antenna operation 


from a systems perspective. 


E. BEAM COUPLING LOSSES VS. NUMBER OF BEAMS 

In the preceding sections, we have examined the effect of various antenna design 
parameters such as sidelobe level and beamwidth on beam coupling losses. At this 
juncture it is important to also investigate the effect of fixing design parameters and 
determining the relationship between beam coupling losses and the number of beams. For 


this investigation, only three scenarios will be investigated 


1. Uniform Weighting 
2. Cosine on a 10 dB Pedestal Weighting 


3. Taylor 25 dB, n = 6 Weighting 
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For each one of the amplitude weights listed, the beam coupling losses for each crossover 
level were compared for the center beams. Figure 14 shows a graph of the coupling loss 
versus the number of beams for a 25 dB Taylor distribution with a crossover level of 3 
dB. The difference in coupling loss between the 2 beam case and the 6 beam case is less 
than 0.9 dB. The results indicate that the addition of a beam has little or no effect if its 
scan direction is greater than several beamwidths from the center beam. When this 
condition is true, the new entries in the 8 matrix due to the added beam are small, and 


the eigenvalues are primarily determined by the elements near the diagonal. 
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V. CONCLUSIONS 


Analysis of the data presented in this report supports several conclusions. These 


conclusions are 


1. For arrays, the pertinent loss equations are readily implemented and solved using 
commercial software packages. Availability of these packages makes it possible 
for personnel involved with array antennas (especially junior officers) to more 
completely understand why the equipment behaves as it does. A more detailed 
understanding of equipment performance translates to increased mission readiness 
and improved program management. 


2. Beam coupling losses are dependent upon the basic array antenna design 
parameters and therefore must be traded off with beamwidth, sidelobe level, and 
gain requirements. 


3. The reduction in array gain observed when using multiple beam techniques varies 
as a function of amplitude taper and crossover level, and is related to radiation 
efficiency. The reduction is not overwhelming, in most cases, neither is it 
insignificant. 


4. The curves of the coupling loss vs crossover were distinctly different for low 
sidelobe and uniform arrays. 


5. Future research should examine planar arrays and include the effect of the 
radiating element factor. The option of integrating the radiation patterns 
numerically (rather than summing the aperture coefficients) should also be 
considered. 

Array antennas are in use throughout the radio frequency portion of the spectrum, 


and recent printed circuit arrays operate at frequencies as high as 800 Gigahertz (GHz) 


[Ref. 10]. Phased array technology will continue to grow in use due to its 
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advantages in sidelobe suppression, scanning speed, and compactness. It is necessary to 
understand the benefits and limitations of this technology if it is to be exploited to its full 


potential. 
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APPENDIX A - CALCULATED DATA 
A. 2 BEAM DATA 
TABLE 2. 2 BEAM UNIFORM - 20 ELEMENTS 
Crossover Level (dB) q’ (dB) 


6.6000 6.0099 0.5000 3.0103 


a 
| 
aaa 


ed ii 
= oe 
| a 
= eter tia 


13.4200 0.2168 0.8218 0.8522 





at 


TABLE 3. 2 BEAM UNIFORM - 100 ELEMENTS 


a a a 
eo 
ithe 
a a al 
ht 
Se ea a 
Sto 
aes 













2.6650 | Oma] : 0.8232 0.8448 
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TABLE 4. 2 BEAM COSINE (10 dB PEDESTAL) - 20 ELEMENTS 


7.8200 0.3289 0.6716 
9.5200 0.3307 0.7515 
1.8900 0.2182 0.8209 a 


oe 
= 
a 
ne 


16.3500 | 0.0152 0.9850 : 0.0655 
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TABLE 5. 2 BEAM COSINE (10 dB PEDESTAL) - 100 ELEMENTS 


SS 

A: 1.0000 0.8909 
a 
ee 
oe 
Se pt 





q’ (dB) 
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TABLE 6. 2 BEAM COSINE (15 dB PEDESTAL) - 20 ELEMENTS 


= = 
a 
= 
= = 
a 


17.5900 0.0030 | 0.9970 0.0130 
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TABLE 7. 2 BEAM COSINE (15 dB PEDESTAL) - 100 ELEMENTS 












Crossover Level (dB) q’ (dB) 


et ed 

1.0981 2.3279 
1.5482 1.7754 
1.8723 1.2858 
2.1425 0.2374 0.9249 
2.3587 0.6165 
Sl esl al Mical 
cl ae a Mal 


10 i® 3.1696 0.0080 0.9920 0.0348 
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TABLE 8. 2 BEAM TAYLOR 20 dB n = 4 - 20 ELEMENTS 


= sn 
SL a 
a 
= 
= a 














13.7900 0.0043 


ne Ben 14.5000 0.0109 0.9892 0.0471 


15.1300 0.0146 0.9856 0.0629 
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TABLE 9. 2 BEAM TAYLOR 20 dB n = 4 - 100 ELEMENTS 


Fe 2.4307 0.0605 0.9429 0.2552 
Si at el Ml Mal 
SL 
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TABLE 10. 2 BEAM TAYLOR 30 dB n = 5 - 20 ELEMENTS 


8.1800 0.5254 0.6556 1.8338 
ae 9.9500 0.377 0.7258 1.3915 
11.4200 0.2689 0.0391 1.0343 


Khia 
= ae 
= a 
Ge 
ro =| 


16.5400 


17.3100 
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TABLE 11. 2 BEAM TAYLOR 30 dB n = 5 - 100 ELEMENTS 


ts 
10 


3.0074 0.0482 0.9540 0.2044 | 





aye. 


TABLE 12. 2 BEAM TAYLOR 40 dB n = 6 - 20 ELEMENTS 


10 19.1700 0.0362 





54 


TABLE 13. 2 BEAM TAYLOR 40 dB n = 6 - 100 ELEMENTS 


Crossover Level iia a a q’ (dB) 


a 
ny Soe li el 


ies =] | 
Gaia Sl Mell al 
a il lel 


3.0074 0.0574 | 0.9458 0. 2422 





2p) 


B. 4 BEAM DATA 


TABLE 14. 4 BEAM UNIFORM - 20 ELEMENTS 


aid 
= hel nla 

[| se] la 
[=| lle 
> = ae 
| te) 
[el 
dh 
MN inal ese 
=e 


26.53 63.31 | 0.06 | 101.08 1.43 10.72 
















Notes: 
L. = a 
6.5 4 B» 5 ow B34 


a Ad,., = Adg,3 = Ads, 
os Bis = Br, 


Ad,.; = Agr, 
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TABLE 15. 4 BEAM UNIFORM - 100 ELEMENTS 


at par pal a 

a a i a 
=== =P = 
2S ee 
ae pla 


i i i al 
ll 
he fh a 
is is Me tl 


Orch? al i 07054 "Or 2 79 OUD) lS: 525 | eer 0: 722 | 


Notes: 





i Bis = Box = B34 
_ Ad,., = Ag,3 = Adg3, 
Bis = Br 


3 Ad... = Adz; 
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TABLE 16. 4 BEAM COSINE (10 dB PEDESTAL) - 20 ELEMENTS 


ata l 
== 
a 
te 
= tener 
tr 
ol a 
a 
Se es 


Notes: 











ies ov B.; = B34 
* Ady.» =Ad,; = Ads 
“oe = Br4 


*Ady.3 =Ao4 
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TABLE 17. 4 BEAM COSINE (10 dB PEDESTAL) - 100 ELEMENTS 





ce ae 
| 
papa pat pe 
| i i a a 
Pe i a Ka 
| Yi es Mc il 
ee 
i i i i is ai 
: 3 










5.456 11.673} 0.009 17.667] 0.04 









18.576] 0.08 |0.982 


0.024 6.349 0.012 — Oe all ag 


Notes: 





Bi. = Boz = Bs4 
; Adg,.. = 4¢23 = Ag, 
4 Bi3 = Bos 


f Ad,.3 = Ado, 


59 


TABLE 18. 4 BEAM COSINE (15 dB PEDESTAL) - 20 ELEMENTS 


0:28 






0.01 


Notes: 
Bia = Bos = Bs 
* Agia = Adr3 = Ads, 
* Bis = Bra 


i Adg,.3 = Ady, 
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TABLE 19. 4 BEAM COSINE (15 dB PEDESTAL) - 100 ELEMENTS 


ial Go 
tea ie 
el a 

Soe a 
eee |= 
so is i iad i 


0.024 | 6.168 | 0.003 |12.409 | 0.002 | 18.804 0.960 































——-.=- | 


13.073 | 0.001 


= — ae —_— —_———— a 1 = ie <—- 


10.003 | 6.494 | -0.002 














0.996 | 





19.834. 





0.02 


eee == 
Notes: 
" Bia = Bos = Bas 
* Ady. = Adrs = Ads 
* Bis = Brg 
* Ady; = Ad), 
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TABLE 20. 4 BEAM TAYLOR 20 dB n = 4 - 20 ELEMENTS 


Riad eM ed ea a Dd 
pt | | Se 
[se 
ps fae ae 
MEM RDEC 
a 
a 
a 
a 
a 

en 

Breas area 


10 91. 59 0.00 |107.30 .27 10.94 







oN 
© 
WS) 
nm 
S 
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a NR am 






Notes: 
' Biz = Bos = Bas 
* Adi, = Ady; = Adss 
* Bis = Bos 


" Ad,; = Ado, 
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TABLE 21. 4 BEAM TAYLOR 20 dB n = 4 - 100 ELEMENTS 


Ce KN A cd 
ee periort tse] =| | 
eater) P=) =| =| - 
Soe 
ee Peale ol || =| 
eee Peale] f=] =| = 
So Gana ame 
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0.017 5.480 | -0.028 | 11.012 0.038 | 16.650 








ee ~~ = 





———— == 











1P559 le “O22 |>r0e95 





0.051 0.022 





5.770 





0.008 11.600 





Notes: 
' Biy = Boz = Bs4 
* Adi. = Adr3 = Adsy 
* Bi = Bo4 
* Adi; = Ady, 
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TABLE 22. 4 BEAM TAYLOR 30 dB n = 5 - 20 ELEMENTS 


ae 
a 
= 
= = 
= 
= 
a 
| =| |e ae 
=| 
7 














0.04 34.70 ). 94.81 -0.0 110.83 


Notes: 
' Bir = Boz = Bsa 
‘ Ag. = Ag,; = Ad, 
* Bis = Boy 
* Adi; = Ady, 
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TABLE 23. 4 BEAM TAYLOR 30 dB n = 5 - 100 ELEMENTS 


eps Te le [ale 
See Le pap [= 
a a ae 
2G See 
ae eee ae | =| 
i lal ll 
a ae 
Gs 














0.006} 19.375] 0.54] 0.884 












0.40; 0.911 





0.007 


20.409 





Notes: 


: Bir = Br3 = Bs. 
. Ad,.. = Ag; = Adz, 
. Bis = Bo4 


: Ad,.3; = Ad 4 
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TABLE 24. 4 BEAM TAYLOR 40 dB n = 6 - 20 ELEMENTS 


1 0.75 12.78 0.29 26.25 0.03 41.57 3.98 


3 0.41 2 -~0.01 2ye2 
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TABLE 25. 4 BEAM TAYLOR 40 dB n = 6 - 100 ELEMENTS 


al I 
JSS Ce 
JSS SeSee 
ete reel mien| = 
A Gls Ss a Se 
ieee ee) exe=| [=| P= 
SSS eae 
ES ae 
Sooo 
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Notes: 

















‘ Bi. = Br3 = Bs4 
* Ady. = Ada; = Adz, 
4 Bis = Br, 


. Ad,.. = Ad? 
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APPENDIX B. 4 BEAM RADIATION EFFICIENCY 


This Appendix describes the radiation efficiency calculation for the case of 4 


identical beams. For 4 beams, the determinant for the radiation efficiency factor becomes 


I-h) Bio Bis Big 
Bs, (1-A) B53 Bo, . (49) 
B5, B5, (1-1) B5, 

| By, By, By3 (1-A)| 


Recalling the symmetry of the coupling coefficients, this determinant can be written as 


1-2) By By Ba 
Ben (1-A) By oe 0 (50) 


B,, ue (1-A) B43 
Bar Bio B43 (1-1) 


The evaluation of the determinant analytically proceeds by cofactor expansion. Expanding 


and collecting terms leads to the final result of 


44-443 +(6-K)A? +2(K-A-2)A +(2A+D+1) = 0. (51) 


The following definitions are made to simplify the expression: 
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a=6,,B,, Bs) Y =B,.B4, By 

C= Be nboaibai €=B4,B 5, B41 (52) 
M =Bs)P42Pa2Pa1 %=8o) By Bar B43 

p = Ba, ie ye Bas 


K = |By5/?+ [Bag + [Bal + [Bail + Bail? + [Ba (53) 
C= 1B, 17 1B,;|° + 1B al" 1B si l° + 1B 39171 B41 7 


A = Re{a]+Ref{y]+Re[6] + Refe] 
B = Re[n] + Re[p] + Refo] (54) 
D=C-2B 


Since the coefficients of the fourth order polynomial obtained in (33) are real, the 
Fundamental Theorem of Algebra guarantees exactly four complex roots. Additionally, 
the roots of the polynomial correspond to the eigenvalues of the beam coupling matrix. 
Since the matrix is Hermitian, all eigenvalues (and therefore the polynomial roots) are 
guaranteed to be real. The largest root (dominant eigenvalue) is substituted into equation 


(38) to obtain the maximum radiation efficiency for the system. 
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APPENDIX C. COMPUTER SOURCE CODE 


A. INTRODUCTION 

The MATLAB™ (Matrix Laboratory) software package by The MathWorks® offers 
efficient implementation of the equations describing the coupling interactions and 
radiation efficiency. Since it is a programming language based on vector/matrix 
manipulation, it is fairly straightforward to enter the desired array coefficients, and have 
the computer calculate the solution. Additionally, many built in routines provide 
flexibility without the necessity of laboring to produce FORTRAN or PASCAL code. 
These features allow the programmer to concentrate primarily on the behavior and 
performance of the algorithm. The matrix nature of MATLAB also allows more efficient 
computation when data are presented in vector format. Acceptable throughput rates were 
obtained on a 386SX personal computer with an installed math co-processor. 

For this thesis, the computer coding proceeded in a modular fashion. Principal 
tasks were coded as functions, with the main procedure calling the functions through 
arguments. Functions written included 

1. Cosine on a Pedestal amplitude weight generator. 
2. Taylor amplitude weight generator. 
3. Beam coupling loss generator. 


4. Array Pattern function 
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The source code for these functions is included in this Appendix. These functions were 
the principal calculational tools utilized in computing the raw data. Sufficient 
documentation is included within the function M-Files to allow straightforward analysis 


of program flow. 


B. MATLAB SOURCE CODE 


1. COSINE ON A PEDESTAL AMPLITUDE WEIGHTS 


function y =cospedwn(N,p) 

% function Y=COSPEDWN(N,p) 

% 

% 

% This function calculates amplitude weights for a cosine on a pedestal 

% window. The input arguments is the number of array elements, (N) and 
% the pedestal size in dB. The function returns the amplitude coefficients. 
% 

% 

clg; 

ped = 10.*(-p/20) 

1=1:N; 

xl=2*(1-1)/(N-1) -1; 

y = (((l-ped)*abs(cos(xl*pi/2))) + ped).’; 


2. TAYLOR AMPLITUDE WEIGHTS 


function a=taylor(slr,n_bar,N) 

% 

% function a=TAYLOR(sir,n_bar,N) 

% 

% This function calculates amplitude weight coefficients for an 

% n element array based on a Taylor amplitude distnbution. The 

% input arguments are the desired main lobe - side lobe ratio (dB), 
% the sidelobe parameter n_bar, and the number of array elements. 
% 

% The output variable is the amplitude weights an. 

% The procedure is outlined in the text " The Handbook of Antenna 
% Design (vol 2)” edited by A.W. Rudge, K. Milne, A.D. Olver, and 
% P. Knight. 

% 

clg; 

p=n_bar-l; 

dbamp =20/log(10); 

sll=exp(abs(slr)/dbamp); % Calculate sidelobe level as a ratio 
A=(1/pi)*log(sll+sqrt(sll.~2-1)); % Generate intermediate parameters 
sigma =n_bar/(sqrt(A.-2 + (n_bar-1/2).*2)); 

f=ones(p,!); 


% Begin calculation of amplitude weights 
for m=1:p 


7A 


z(m) = sigma*sqrt(A.~2 +(m-1/2).*2); 
K(m)=(fact(p)).°2 ./ (fact(p +m) .* fact(p-m)); 


end 


for n=1:p 

P(n) =prod(1-(n.*2 ./ z.~2)); 
end 
f=K.*P; 
J=zeros(N/2,p); 
for 1=1:N/2 

for k=I:p 

J(i,k) = f(:,k). *cos(k. *p1.*2*(i-1)/N); 

end 
end 
g=1+2*sum(J.’); 
% Only 1/2 of the coefficients need to be calculated due to symmetry 
a=[fliplr(g) g].’; 
a=a/max(a); % Normalize the amplitude to max of | 


3. BEAM COUPLING LOSS CALCULATOR 


function beta =couploss(N,crosslvl,win) 

% 

% function beta=COUPLOSSV(N, crosslvl,win) 

% This function calculates coupling losses for overlapping beams 

% in a multiple beam antenna. The multiple beams are generated using 

% a linear array with various choices for amplitude windows. 

% A switch parameter is used to choose the amplitude weighting as follows 


% win == Rectangular window 

%o win == 2 _ Cosine on a pedestal (Pedestal size is input) 
% win == Taylor (Sidelobe level and n_bar are input) 
% 

clg; 


% Initialize parameters 
dB=crosslv]; 

k=2*pi; 

d=.5; 

kd=k*d; 

zeta=(N + 1)/2; 


% Get amplitude coefficients 
if win == 
al =boxcar(N); 
a2=al.’; 
elseif win == 
p=.3199; 
al =cospedwn(N,p); 
a2=al.’; 
elseif win == 
slr=20; 
n_bar=4; 
al =taylor(slr,n_bar,N); 
a2=al.’; 


end 
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% Get crossover squint angles 
level=10 .*(-dB/20); 

u=Q; 

xi=1; 

delta_u=1/10000; 
a=al(N/2+1:N).’; 

m= 1:N/2; 

Max =a*ones(10, 1) +a*ones(10, 1); 


while xi > = level 
u=ut+delta_u; 
z1 =a*exp(j*(2*m-1)/2*kd*u).’; 
z2=a*exp(-j*(2*m-1)/2*kd*u).’; 
xi= 1/Max*(z1 +22); 
end 
u_squint=u; 
crossover _des=level; 
crossover_comp = xi; 


% Initialize scan angles and sums 
up =kd*u_ squint; 

uq =-kd*u_squint; 

du=uq-up; 

dtheta = (acos(du/kd))* 180/pi; 


% Set up matrices and calculate summations 


Al =al*a?2; 
[=zeros(N,N); 
a= 72—0:73=—0; 
for m=1:N 
for n=1:N 
if n==m 
I(m,n) =2*kd; 


Z1=Z1+Al1(m,n).*1(m,n).*exp(-)*m*du); 
Z2=Z2 + Al(m,n).*1(m,n). *exp(j*0); 
Z3=Z3 + Al(m,n).*1(m,n)*exp(j*0); 
else 
I(m,n) =2*((sin((n-m)*kd))/((n-m))); 
Z1=Z1+A1(m,n).*1(m,n).*expG*(up*m-uq*n)); 
Z2=Z2+Al1(m,n).*1l(m,n).*expg*up*(m-n)); 
Z3=Z3+Al1(m,n).*1(m,n).*expG*uq*(m-n)); 
end 
end 
end 
% Calculate coupling coefficient here 
ex=1.*exp(j*zeta*du); 
X=ex.*Z1; 
be=((Z2). (.5)) .* ((Z3).°(.5)); 
b=xX./Y; 
beta =abs(b); 


4. ARRAY PATTERN FUNCTION 


function x=arrpat(a,psi_primel,n) 
% 

% X=ARRPAT(a,psi_primel,n) 
% 


% This function calculates the far-field array pattern for a 
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% linear array. Input arguments are the amplitude weight coefficients 
% the desired steering angle off broadside, and the number of elements. 
% The interelement spacing is automatically assumed to be lambda/2. 
% 
clg; 
if psi_primel <= 90 % Refer steering angle to broadside. 
psil =(pi/2) - psi_prime]l *(pi/180); 
else 
psil=psi_primel *(pi/180), 
end 
[j,k] =size(a); 
ifj > 
a=a.’; 
end 
lambda = 1; % Normalize wavelength 
k = 2*pi/lambda; 
d=lambda/2; % Interelement spacing 
kd=k*d; 
phil = -kd*cos(psil); % Generate phase shift terms to steer the 
% beam. 
if rem(n,2) ~ = 0 
z=-(n-1)/2:(n-1)/2; 
else 
z= linspace(-n/2,n/2,n); 
end 
psi=linspace(0,2*pi,256); 
el =exp(j*z*phil); % Phase weighting terms 
se=size(el); 
% Generate amplitude and phase weighted terms 
al =a(n/2+1:n); 
X1=zeros(256,n/2); 
for 1=1:256 
for z = 1:n/2 
X1(1,2) =(al(:,z).*cos(kd.*cos(psi(l)).*(2*z-1)/2 + ((2*z-1)/2.*phil))); 
end 
end 
ml =2*abs(sum(X1.’)); % Sum the elements to get the array factor 
x=ml1/max(ml1); 
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